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ABSTRACT 

This paper presents Spitzer-IRS A/AA ^ 600 spectroscopy of the CO2 15.2 ^m bending mode toward 
a sample of 50 embedded low-mass stars in nearby star- forming clouds, taken mostly from the "Cores 
to Disks (c2d)" Legacy program. The average abundance of solid CO2 relative to water in low-mass 
protostellar envelopes is 0.32 ±0.02, significantly higher than that found in quiescent molecular clouds 
and in massive star forming regions. It is found that a decomposition of all the observed CO2 bending 
mode profiles requires a minimum of five unique components. In general, roughly 2/3 of the CO2 ice 
is found in a water-rich environment, while most of the remaining 1/3 is found in a CO environment 
with strongly varying relative concentrations of CO2 to CO along each line of sight. Ground-based 
observations of solid CO toward a large subset of the c2d sample are used to further constrain the 
C02:C0 component and suggest a model in which low-density clouds form the C02:Il20 component 
and higher density clouds form the C02:C0 ice during and after the freeze-out of gas-phase CO. The 
abundance of the C02:C0 component is consistent with cosmic ray processing of the CO-rich part 
of the ice mantles, although a more quiescent formation mechanism is not ruled out. It is suggested 
that the subsequent evolution of the CO2 and CO profiles toward low-mass protostars, in particular 
the appearance of the sphtting of the CO2 bending mode due to pure, crystalline CO2, is first caused 
by distillation of the C02:C0 component through evaporation of CO due to thermal processing to 
~ 20 — 30 K in the inner regions of infalling envelopes. The formation of pure CO2 via segregation 
from the II2O rich mantle may contribute to the band splitting at higher levels of thermal processing 
(> 50 K), but is harder to reconcile with the physical structure of protostellar envelopes around 
low-luminosity objects. 
Subject headings: astrochemistry — circumstellar matter — dust, extinction — ISM: evolution 



1. INTRODUCTION 

Although CO2 is not an abundant gas-phase molecule 
in molecular clouds, it is one of a small number of molec- 
ular species consistently found in very high abundances 
inside ice mantles on dust grains (> 10~^ with r espect 
to H9 iGerakines et all [l999l : IWhittet et al.ll2007D . The 
generally high abundance of solid CO2 became apparent 
with the spectroscopic surveys c onducted with the In- 
frared Space Observatory (ISO) (jGerakines et al.lll999t 
iNummclin et al. 2001). Other species known to belong 
to this class of very abundant molecules are CO and 
H2O. In less than 10% of dark cloud regions surveyed, 
methan ol (CII3OH) is also found with similar abun- 
dances (i Pontoppidan et al.l[2003at iBoogert et al.ll2007f ). 
Depending on the density and temperature of a cloud, 
the CO is found partly in the gas-phase and partly frozen 
onto grain surfaces, while t he CO2 and II2O a re com- 
pletely frozen as ice mantles (jBergin et al.lll995l ). except 
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in very hot or shocked regions ("Boonm an et al.l l2003t 
iNomura fc Milla r 2004; Lahuis et al. 20011). The system 
of CO, CO2, II2O and, under some conditions, CII3OII 
therefore represents the bulk of solid state volatiles in 
dense star forming clouds, and interactions between these 
four species can be expected to account for most of the 
solid state observables. Other species with abundances 
of less than 5% relative to water, such as CII4, NII3, 
OCN-, HCOOH and OCS wiU be good tracers of chem- 
istry and their local molecular environment, but are un- 
likely to strongly affect the molecular environments, and 
therefore the band profiles, of the four major species. 

The formation mechanism of solid CO2 in the cold 
interstellar medium is still not understood, although a 
number of plausible scenarios have been proposed. Since 
the direct surface route, CO -f O ^ CO2, is thought 
to possess a significant activation barrier, it was ini- 
tially suggested that strong UV irr adiation was needed 
to pr oduce the observed CO2 ice (jd'Hendecourt et ahl 
Il985| ) . Laboratory simulations of interstellar ice mixtures 
of H2O and CO confirmed that C O2 is indeed readily 
formed during strong UV photolysis (jd'Hendecourt et al.l 
Il986[ ). and initial detection of abundant CO2 ice around 
UV-luminous massive young stars seemed to confirm this 
picture. However, recent detections of similar abun- 
dances of CO2 in dark clouds observed along lines 
of sight toward background stars, f a.r away fro r n any 
ionizi ng source (iBergin et al.l l2005t iKnez et al] l2005l : 
Whittet et all l2007f) . argue against a UV irradiation 
route to CO2, at least through enhanced UV from nearby 



protostars. iPontoppidanI ()2006[ ) showed evidence for an 
increasing abundance of CO2 ice with gas density in at 
least one low-mass star forming core. Furthermore, the 
original p remise of a barrier to oxygenation of CO is now 
in doubt (jRoser et al.ll200ll ). Consequently, both theo- 
retical and laboratory efforts to understand the forma- 
tion of CO2 are still very active. 

Extensive surveys of the 3.1 and 4.67//m stretch- 
ing mode of H2O and CO ices have been car- 
ried out in a range of differe nt star-forming en- 
vironments (IWhittet et all 119881 : IChiar et all Il995t 
iPontoppidan et al.ll2003b( l. However. CO2 can only be 
observed from space, and surveys have until recently 
been limited to a small sam ple of very luminous young 
stars (iGerakines et al.l ll999V While the 4.27 ^m stretch- 
ing mode of CO2 was detected toward a few background 
field stars by ISO (jWhittet et al.lll998t iNummelin et all 
l2001h . recent Spitzer observations of the 15.2 /zm bend- 
ing mode of CO2 have extended the sa mple of CO2 along 
quiescent liri es of sight cons iderably (jKnez et al.l 120051 : 
iBerginet al]|2005l : IWhittet et al. 20m . 

The profiles of the CO2 ice bands observed in quies- 
cent regions and along lines of sight toward luminous 
protostars show intriguing differences. In particular, 
ice in massive star-forming regions, which presumably 
traces more processed material, tends to show a split- 
ting of the 15.2 ^m bending mode. This splitting has 
been identified as a general property of crystalline pure 
CO2 and is readily reproduc ed in laboratory s i mula- 
tions of in terstellar ices (e.g. lEhrenfreund et al] Il997l : 
Ivan Broekh uizcn et al. 2006). It seems implausible that 
this pure CO2 layer would form directly through gas- 
phase deposition and subsequent surface reactions, and 
it has been suggested that the CO2 segregates from a 
C02:H20: CH30H=1:1 :1 mixture upon strong heating 
(jGcrakinc s et al.l Il999l ) . Anne aling in a laboratory sct - 
ting produces the same effect (jEhrenfreund et anil997f ). 
and strong heating is not unreasonable in the envelopes of 
massive young stars with luminosities in excess of 10^ Lq. 

In this paper, a survey of the 15.2 /xm CO2 bend- 
ing mode toward ^50 young low mass stars using the 
high resolution mode of the Infrared Spectrograph on 
board the Spitzer Space Telescope (Spitzer-IRS) is pre- 
sented. The sample stars have typical luminosities in the 
range 0.1 — 10 Lq, thus bridging the observational gap be- 
tween the background stars and the massive protostars 
from the ISO sample. The Spitzer spectra are comple- 
mented by ground-based observations of II2O and CO 
ices, where available, as well as the archi val spectra from 
the Infrared Space Observatory used in iGerakines et al.l 
11999). While this paper concentr ates on the region 
around 15 /im, iBoogert et al.l ((2007) discusses the ices 
causing the 5-8 /im absorption complex (henceforth re- 
ferred to as Paper I). The 7.7 /im CH4 and 9.0 /im NII3 

bands are described in two separate papers (Oberg et 
al. 2008, submitted and Bottinelli et al., in prep, respec- 
tively) . 

The central questions that are addressed using the new 
Spitzer spectra of the CO2 bending mode are: 

• What are the differences, if any, between CO2 ice in 
massive and low-mass star-forming environments? 

• What is the average abundance of CO2 in low-mass 



protostellar envelopes compared to lower density 
quiescent clouds? 

• In which molecular environments can solid CO2 be 
found, and what are their relative abundances? 

• Which process forms CO2 in CO-dominated envi- 
ronments when the CO accretes from the gas-phase 
at high densities? 

• How does the component of pure CO2 (as measured 
by the well-known splitting of the bending mode) 
form in low-mass protostellar envelopes? 

• What is the evolution of the CO2 ice and what does 
it tell us about protostellar evolution? 

2. THE INFRARED BANDS OF SOLID CO2 

The infrared vibrational modes of CO2 are known to 
be very sensitive to the molecular environment. Observa- 
tions of the band profiles can determine whether the CO2 
molecules are embedded with water, CO and other CO2 
molecules. Solid CO2 has two strong vibrational modes; 
the asymmetric stretching mode centered on 4.27 /im, 
and the bending mode at 15.2 /im. The stretching mode 
is so strong that it is typically saturated along lines of 
sight through protostellar envelopes, and its ^^C02 coun- 
terpart at 4.38 /im is often used for profile analysis in- 
stead. However, neither stretching modes are covered by 
the spectral range of Spitzer-IRS. 

Fortunately, the CO2 bending mode is an excellent di- 
agnostic of molecular environments. For instance, pure 
CO2 will typically produce a split band in the bending 
mode due to Davydov splitting - a long range interaction 
in crystalline materials. Conversely, CO2 embedded in a 
hydrogen-bonding matrix will produce a broad, smooth 
profile. While these differences are relatively well known, 
it is noted that the different formation scenarios from 
gas-grain chemical models make distinct predictions for 
the molecular environment of the CO2. Thus models fa- 
voring a formation route via OH predict that the CO2 
will be found in a water-dominated matrix. 

Indeed, observations of the CO2 stretching and 
bending modes towar d young massive stars with ISO 
(jGerakines et al.iri999t ) have shown that the CO2 ice is 
dominated by a band consistent with CO2 in a hydrogen- 
bonding environment. They also found that a double- 
peaked component, consistent with a relatively pure, 
crystalline CO2 is generally present at a lower level to- 
ward their sample of massive stars. 

3. OBSERVATIONS 

The spectra of the 15.2 /im CO2 bending mode have 
been obtained using the short-high (SH) module of 
Spitzer-IRS, with a spectral resolving power of A/AA ~ 
600, covering 10-19. 5/im, corresponding to lcm~^. The 
Spitzer spectra were obtained as part of the "Cores to 
Disks" Legacy program (PID 172,179) as well as a few 
archival spectra observed as part of the GTO programs 
(PID 2). All SH spectra from the c2d database that show 
clear detections of the CO2 bending mode have been in- 
cluded. The spectra have been reduced using the c2d 
pipeline from basic calibrated data (BCD) products ver- 
sion S13.0.2. For each spectrum, clearly deviant points 
were removed and individual orders were scaled by small 



factors to align the overlapping regions between orders. 
The overlapping regions between these two orders usu- 
ally match very well, which lends support to the reality 
of small features in the spectra for most sources. A small 
number of sources show evidence for absorption by the Q- 
branch of gas-phase CO2 at 15.0 micron. Since this sur- 
vey is concerned with solid CO2, pixels affected by gas- 
phase absorption have been removed from the fits of IRS 
46, WL 12 and DG Tau B. Some of the most embedded 
sources have saturated or nearly saturated bending mode 
bands. For these sources, special care has to be taken to 
ensure that the background level is well-determined. 

In addition to the Spitzer data, the five highest-quality 
spectra of massive YS Os observed with the ISO-SWS 
(jGerakines et all 11999^ have been included. The ISO 
spectra provide a useful comparison of the structure of 
CO2 ices in the warmer, more energetic envelopes of 
young massive stars to the comparatively quiescent en- 
velopes of low-mass YSOs. 

Finally, to study the relation of the CO2 ices with CO 
and water, ground-based spectra of the 4.67 /im stretch- 
ing mode of solid CO and the 3.08 /im stretching mode 
of water ice have been collected using the Infrared Spec- 
trometer and Array Camera (ISAAC) on the Very Large 
Telescope (VLT)^ Mo st of the 4.67Mni VLT-ISAAC 
spectra are pubUshed in lPontoppidan et al.l ([2003b), but 
a significant fraction are previously unpublished spectra 
obtained with NIRSPEC at the Keck Telescope. Sources 
with no CO ice data are typically too faint for useful 
ground-based 4.67 /xm spectroscopy. The column densi- 
ties of water ice are derived using the 3.08 /im ground- 
based spectra or taken from Paper I. The data set is 
summarized in Tabled] The new observations of the CO 
ice bands not covered in lPontoppidan et al.l (|2003bD are 
summarized in Table [51 

4. PROFILE DECOMPOSITION 

4.1. Continuum determination 

To directly compare with dust models, the spectra of 
ice absorption bands have to be converted to an opti- 
cal depth scale. This requires that an appropriate con- 
tinuum be defined, a process somewhat complicated for 
CO2 by the location of the bending mode on the blue side 
of the broad silicate bending mode and on the red side 
of the II2O libration band. Unfortunately, not knowing 
the shape of the underlying continuum, this is a prob- 
lem with no unique solution. In this work, continua for 
each spectrum are constructed by fitting a third-order 
polynomial to the spectral ranges: 13 /im ~ 14.7 /im and 
18.2 /im - 19.5 /im. The shape of the blue wing of the 
silicate bending mode is simulated by a Gaussian in fre- 
quency space with center at 608 cm"-'^ and a FWHM of 
73cm~^. The aim is to construct a shape of the con- 
tinuunr that has a negative second derivative under the 
CO2 band. The same procedure has been used for all the 
Spitzer spectra, and the resulting continua are shown for 
each spectrum in Figures |4] through 1161 

4.2. Laboratory data 

A number of laboratory spectra of CO2 ices have been 
taken from the literature. Ice inventories of envelopes 

''' Based on observations made with ESO Telescopes at the 
Paranal Observatory under programme ID 164.C-0605 



around young low-mass have shown that the ices are 
dominated by II2O, CO and CO2, so this study concen- 
trates on systems involving these three species. In some 
regions of low-mass star formation, CH3OII is found in 
large amounts (up to 25% relative to II2O), but this 
seems to affect only a small subset of the sample pre- 
sented here. 

The available laboratory spectra are divided into 
C02:Il20 mixtures and C02:C0 mixtures, each set with 
distinct characteristics. As in the case of all solid state 
features due to abundant molecules, the band shapes can 
be strongly modified by surface modes, depending on 
the s hape distribution of the dust grains (jTielens et al.l 
Il991h . Astronomical spectra can therefore not be directly 
compared to absorbance laboratory spectra. Rather, 
complex refractive indices must be combined with a dust 
model to calculate opacities relevant for the small irreg- 
ular dust grains of the interstellar medium. This study 
is consequently restricted to laboratory experiments for 
which optical constants have bee n calculated. 

Ehrcnf reund et a l. (1997| ); iDartois et all (119991 ) 
present optical constants for a wide range of C02:C0 
mixtures, as well as a few C02:H20 mixtures, obtained 
under high vacuum (HV) conditions and 2.0 cm~^ 
resolution. More recently, a number of detailed 
studies of relevant C02-rich ices were performed by 
Ivan Broekhuizen et al.l (|2006D and lOberg et all (|2007af) . 
also under HV conditions. While these studies do not 
provide optical constants directly, they do report ab- 
sorbance spectra as well as approximate ice thicknesses, 
making it possible to derive optical constants using the 
Kramers-Kronig relations. 

For the Il2 0-rich ices, the CO2:H2O =14:100 mixture 
at 10 K from lEhrenfreund et al.l (jl997f ) was chosen. The 
CO-rich ices show a CO2 bending model profile that is 
dependent on the mixing ratio. A function is therefore 
constructed that returns an ice spectrum for any relative 
concentration between C02:CO=l:4 and C02:C0=1:1 
by interpolating between the available 10 K laboratory 
spectra within this range. At very low concentrations of 
CO2 relative to CO (CO2:CO<l:10), the bending mode 
becomes quite narrow, and the shape becomes indepen- 
dent of concentration. At very high concentrations of 
CO2, the bending mode exhibits the well-known split, 
characteristic of a crystalline structure of the ice. The 
peaks are very narrow - of order 1 cm~^, and so the hi gher 
resolution data from Ivan Broekhuizen et al.l ()2006f ) are 
used for pure CO2 ice. This spectrum, however, suf- 
fers from a misalignment in the spectrometer optics at 
0.5 cm~^ resolution so that the bending mode is too weak 
by a factor of 3 relative to the stretching mode and 
the noise is relatively high. As a result, it was nec- 
essary to scale the absorbance of the bending mode to 
fit with the band st r ength of 1.1 x lO^^^cm^^ reported 
by iGerakines et al.l (|1995f ) before calculating the corre- 
sponding set of optical constants. The noise in the spec- 
trum is reduced by fitting a number of Gaussians, rather 
than smoothing it, which would reduce the resolution. 
Because all the laboratory spectra used were obtained 
under high-vacuum conditions, they may be contami- 
nated by II2O. 

At 15 /im the dust grains are most likely well into the 
Rayleigh limit of 27ra <C 15 /im, where a is the radius of 



TABLE 1 

Observed sample of embedded young stars 



Source 



CO2 



CO and H2O RA [J2000] DEC [J2000] Observation ID 



W3 IRS5 


ISO-SWS 


NIRSPEC 


02 25 40.8 


+62 05 52.8 


42701302 


L1448 IRSl 


Spitzer 


NIRSPEC 


03 25 09.4 


+30 46 21.7 


5656832 


L1448 NA 


Spitzer 


- 


03 25 36.5 


+30 45 21.4 


5828096 


L1455 SMMl 


Spitzer 


- 


03 27 43.2 


+30 12 28.8 


15917056 


RNO 15 


Spitzer 


NIRSPEC 


03 27 47.7 


+30 12 04.3 


5633280 


IRAS 03254+3050 


Spitzer 


NIRSPEC 


03 28 34.5 


+31 00 51.2 


13460480 


IRAS 03271+3013 


Spitzer 


NIRSPEC 


03 30 15.2 


+30 23 48.8 


5634304 


Bl a 


Spitzer 


NIRSPEC 


03 33 16.7 


+31 07 55.1 


15918080 


Bl c 


Spitzer 


- 


03 33 17.9 


+31 09 31.0 


15916544 


IRAS 03439+3233 


Spitzer 


NIRSPEC 


03 47 05.4 


+32 43 08.5 


5635072 


IRAS 03445+3242 


Spitzer 


NIRSPEC 


03 47 41.6 


+32 51 43.8 


5635328 


L1489 IRS 


Spitzer 


ISAAC 


04 04 42.6 


+26 18 56.8 


3528960 


DC Tau B 


Spitzer 


NIRSPEC 


04 27 02.7 


+26 05 30.5 


3540992 


GL 989 


ISO-SWS 


NIRSPEC 


06 41 10.2 


+09 29 33.7 


72602619 


HH46 IR 


Spitzer 


ISAAC 


08 25 43.8 


-51 00 35.6 


7130112 


CED 110 IRS4 


Spitzer 


- 


11 06 46.6 


-77 22 32.4 


5639680 


B35 


Spitzer 


- 


11 07 21.5 


-77 22 11.8 


5639680 


CED 110 IRS6 


Spitzer 


ISAAC 


11 07 09.2 


-77 23 04.3 


5639680 


IRAS 12553-7651 


Spitzer 


- 


12 59 06.6 


-77 07 40.0 


9830912 


ISO Chall 54 


Spitzer 


- 


13 00 59.2 


-77 14 02.7 


15735040 


IRAS 13546-3941 


Spitzer 


- 


13 57 38.9 


-39 56 00.2 


5642752 


IRAS 15398-3359 


Spitzer 


- 


15 43 02.3 


-34 09 06.7 


5828864 


CSS 30 IRSl 


Spitzer 


ISAAC 


16 26 21.4 


-24 23 04.1 


5647616 


WL 12 


Spitzer 


ISAAC 


16 26 44.2 


-24 34 48.4 


5647616 


GY 224 


Spitzer 


NIRSPEC 


16 27 11.2 


-24 40 46.7 


9829888 


WL 20 


Spitzer 


- 


16 27 15.7 


-24 38 45.6 


9829888 


IRS 37 


Spitzer 


ISAAC 


16 27 17.6 


-24 28 56.5 


5647616 


IRS 42 


Spitzer 


ISAAC 


16 27 21.5 


-24 41 43.1 


5647616 


WL 6 


Spitzer 


ISAAC 


16 27 21.8 


-24 29 53.3 


5647616 


CRBR 2422.8-3423 


Spitzer 


ISAAC 


16 27 24.6 


-24 41 03.3 


9346048 


IRS 43 


Spitzer 


ISAAC 


16 27 27.0 


-24 40 52.0 


12699648 


IRS 44 


Spitzer 


ISAAC 


16 27 28.1 


-24 39 35.0 


12699648 


Elias 32/IRS 45 


Spitzer 


ISAAC 


16 27 28.4 


-24 27 21.4 


12664320 


IRS 46 


Spitzer 


ISAAC 


16 27 29.4 


-24 39 16.3 


9829888 


VSSG 17/ IRS 47 


Spitzer 


ISAAC 


16 27 30.2 


-24 27 43.4 


5647616 


IRS 51 


Spitzer 


ISAAC 


16 27 39.8 


-24 43 15.1 


9829888 


IRS 63 


Spitzer 


ISAAC 


16 31 35.7 


-24 01 29.5 


9827840 


L1689 IRS5 


Spitzer 


- 


16 31 52.1 


-24 56 15.2 


12664064 


RNO 91 


Spitzer 


ISAAC 


16 34 29.3 


-15 47 01.4 


5650432 


W33 A 


ISO-SWS 


ISO-SWS 


18 14 39.7 


-17 52 02.0 


32900920 


GL 2136 


ISO-SWS 


ISO-SWS 


18 22 27.0 


-13 30 10.0 


33000222 


Serp S68 


Spitzer 


- 


18 29 48.1 


+01 16 42.5 


9828608 


EC 74 


Spitzer 


NIRSPEC 


18 29 55.7 


+01 14 31.6 


9407232 


SVS 4-5 


Spitzer 


ISAAC 


18 29 57.6 


+01 13 00.6 


9407232 


EC 82 


Spitzer 


ISAAC 


18 29 56.9 


+01 14 46.5 


9407232 


EC 90 


Spitzer 


ISAAC 


18 29 57.8 


+01 14 05.9 


9828352 


SVS 4-10 


Spitzer 


ISAAC 


18 29 57.9 


+01 12 51.6 


9407232 


CK4 


Spitzer 


NIRSPEC 


18 29 58.2 


+01 15 21.7 


9407232 


CK 2 


Spitzer 


ISAAC 


18 30 00.6 


+01 15 20.1 


11828224 


RCrA IRS 5 


Spitzer 


ISAAC 


19 01 48.0 


-36 57 21.6 


9835264 


RCrA IRS 7A 


Spitzer 


ISAAC 


19 01 55.3 


-36 57 22.0 


9835008 


RCrA IRS 7B 


Spitzer 


ISAAC 


19 01 56.4 


-36 57 28.0 


9835008 


CrA IRAS 32 


Spitzer 


- 


19 02 58.7 


-37 07 34.5 


9832192 


S140 IRSl 


ISO-SWS 


NIRSPEC 


22 19 18.4 


+63 18 45.0 


22002135 


NGC 7538 IRS 9 


ISO-SWS 


NIRSPEC 


23 14 01.7 


+61 27 20.0 


09801532 


IRAS 23238+7401 


Spitzer 


- 


23 25 46.7 


+74 17 37.2 


9833728 



TABLE 2 
New CO ice band observations 



Source 



r(C0:H20) (red)" '' T(Purc CO) (middle)'^ t(C0:C02) (blue)" 



W3 IRS5 O.llit 0.05 

L1448 IRSl O.lOit 0.05 

RNO 15 0.12± 0.03 

IRAS 03254+3050 0.12± 0.06 



IRAS 03271+3013 


<0.5 


Bl a 


<0.5 


IRAS 03439+3233 


0.09+ 0.05 


IRAS 03445+3242 


0.30+ 0.05 


DG Tau B 


0.13+ 0.05 


GL989 


0.19+ 0.03 


HH46 IR 


0.80+ 0.08 


GY224 


<0.2 


IRS 37 


0.08+ 0.02 


GL 2136 


0.20+ 0.05 


EC 74 


<0.10 


CK4 


<0.03 


S140 IRSl 


0.05+ 0.02 


NGC 7538 IRS 9 


0.52+ 0.05 



0.27+ 0.06 
0.18+ 0.06 
0.27+ 0.03 
0.19+ 0.10 
1.0+ 0.5 
4.0+ 2.0 
0.25+ 0.10 
1.16+ 0.20 
0.12+ 0.10 
0.37+ 0.05 
0.68+ 0.06 
<0.2 

0.45+ 0.04 
0.10+ 0.05 
0.75+ 0.06 
0.48+ 0.05 
0.04+ 0.02 
3.17+ 0.10 



0.06+ 0.05 
0.09+ 0.03 
0.09+ 0.03 
0.11+ 0.05 
<0.5 
<0.3 

0.08+ 0.04 
0.33+ 0.05 
0.11+ 0.05 
0.16+ 0.05 
0.16+ 0.05 
<0.2 

0.10+ 0.05 
<0.05 
0.10+ 0.06 
0.07+ 0.03 
<0.03 
0.30+ 0.03 



" Decomposition as described in lPontoppidan et al.l Il2003bl 'l. 

^ Adopted conversion to column density: AfcO:H20 = 16.0cm~^rco:H2O^~^i where A = 1.110" 
strength. 

'^ Adopted conversion to column density, using a continuous distribution of ellipsoids (CDE): A^co ■ 

"^ AfC0:C02 =3.0cm-lTCO:CO2^"^- 
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Fig. 1. — Com parison of the absorbance spectrum of 
Ivan Broekhuizen et al. : (2006) relevant for pure CO2 ice (red, dash- 
dotted curve) with that calculated for a CDE distribution of par- 
ticles in the Rayleigh limit (blue, solid curve). The dashed curve 
shows the component fit used to eliminate the noise. 

the largest grains. This means that scattering of hght 
out of the hne of sight is unimportant, and the opacities 
can be treated as pure absorption coefficients. A contin- 
uous distribution of eUipsoids (CDE) is used to convert 
the optical constants to opacities. This is a convenient 
method of simulating the effect of irregularly shaped 
grains, found to work well for ice bands in the Rayleigh 
limit and it has been use d successful ly for solid CO 
(iPontoppidan et all l2003bD and CO2 (jGerakines et al.l 
Il999( ). Figure [1] illustrates the process of converting the 
absorbance spectrum of pure CO2 into an opacity that 
can be used for comparing with the Spitzer spectra. The 
laboratory spectra are summarized in Table [3l 

4.3. Component analysis 

The strategy adopted here for analyzing the general 
shape of the CO2 ice bending mode in low-mass young 





TABLE 3 
Laboratory spectra 




Mixture 


TfKl 


Resolution [cm 


^ 


Reference 


CO2:H2O=14:100 


10 


2 






CO2:CO=4:100 


10 


2 






CO2:CO=26:100 


10 


2 






CO2:CO=70:100 


10 


2 






CO2:CO=112:100 


10 


2 






Pure CO2 


15 


0.5 




2 



1 lEhrenfreund et al.l ||T997|') 

2 Ivan Broekhuizen et "all (1200^) 

stellar envelopes is to determine the minimum number of 
unique components required to fit all the observed bands. 
In this context, a unique component is a band that only 
changes its relative depth, but not it s shape from source 
to sour ce. This approach was used in lPontoppidan et al.l 
(|2003b| ) to determine that only three unique components 
could be used to fit the 4.67 /im stretching mode of solid 
CO, and is also used in Paper I to decompose the 5- 
8 /im complex. The three unique CO components were 
1) A broad, red-shifted component associated with CO 
in a water-rich mantle, 2) a component indistinguishable 
from pure CO and 3) a narrow, blue-shifted component 
due to either CO in a CO2 environment, or CO in a 
crystalline form. The three CO components were named 
"red" , "middle" and "blue" , and is seen below, the "red" 
and "blue" components have counterparts in the CO2 
bending modes. Note that the "red" and "middle" com- 
ponents are also often referred to in the literature as "po- 
lar" and "nonpolar" , respectively. The bending mode of 
CO2 probes ice structures that are somewhat more com- 
plicated than those probed by CO. While most of the 
CO ice desorbs efficiently at temperatures higher than 
20 K, the less volatile CO2 ice goes through several ad- 
ditional structural changes. The most characteristic is 
the appearance of the double peak seen in pure CO2 (see 
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Fig. 2. — Empirically determined profiles for the three dominant 
components of the CO2 bending mode. The data are as follows. 
Top: RNO 91 - IRS 42, fit using pure CO2. Middle: CRBR 2422.8- 
3423 - IRS 51, fit using CO2:H2O=14:100. Bottom: IRS 63 - IRS 
42, fit using C02:C0 = 26:100. 

Figure [T]) . 

In order to empirically derive the shape of the compo- 
nents of the CO2 bending mode, pairs of spectra are sub- 
tracted. If the spectra are superpositions of a small num- 
ber of components with varying relative contributions, it 
will be possible to isolate each component. The three 
dominant components determined this way are shown in 
Figure [2] where they are compared to laboratory simula- 
tions. 

Consequently, it is found that the minimum number of 
unique components required to fit all the observed CO2 
bending mode profiles is five: 

• An II2O rich component, modeled with a 
laboratory spectrum with a concentration of 
CO2:H2O=14:100. This is required to fit the red 
wing of almost all the observed bands. This is re- 
ferred to as the "red" component. 

• A component with a roughly equal mixture of CO2 
and CO. Strictly speaking, this is not constructed 
as a unique component since the C02:C0 mixing 
ratio is included as a free parameter. This is re- 
quired because the empirical profiles of the blue 
component, one of which is shown in Figure [21 
have varying widths. This behavior can be repro- 
duced by varying the concentration of CO2 in CO. 
While the band profiles are only available for a set 
of discrete mixing ratios, profiles are constructed 
with arbitrary mixing ratios by linearly interpo- 
lating the available laboratory profiles at each fre- 
quency point. In effect, this allows a measurement 
of the C02:C0 mixing ratio for each observed CO2 
bending mode. This is referred to as the "blue" 
component. 

• A component in which CO2 is very dilute in an oth- 
erwise pure CO ice. This is a very narrow compo- 
nent centered on 15.15/.tm (660 cm~^). In practice. 



this band is modeled by a CO2:CO=4:100 labora- 
tory spectrum. Note that at such high dilution, 
the shape of the band is not sensitive to the exact 
mixing ratio. 

• A component of pure CO2, producing the charac- 
teristic double-peaked structure often seen in pro- 
tostellar sources. 

• An additional, relatively narrow component on the 
red side of the main band. This component is 
unambiguously identified in only a few sources, 
most of them the massive young stars included 
from the ISO sample. The component has in the 
past been identified as an interaction with CII3OII 
in strongly annealed ices with the mixing rati o 
H20:C02:CH30H=1:1:1 (|Cerakines et al.lTl999h . 
Since there is no laboratory spectrum of the shoul- 
der in isolation, the component is modeled empiri- 
cally using a superposition of two Gaussians: 



7'shouldcr(i^)/7'0 — 

-{v — 645 cm~^)^ 



exp 



1.85 exp 



2(2.1cm-i)2 
— {v — 650 cm^^)^ 



2(2.8cm-i)2 
where tq is a scale factor. 



The "red" and "blue" components generally dominate 
the bending mode profiles, and the total CO2 abundance. 
The remaining three components represent subtle differ- 
ences due to trace constituents. It is stressed that all the 
components correspond to distinct and plausible molec- 
ular environments. The relative contributions to three 
typical CO2 bending mode profiles are sketched in Fig- 
ure [S] 

A non-linear least squares fitting routine from the IDL 
library of C. Markwardt * is used to find the best fit 
to each bending mode spectrum. Since three of the five 
components used to interpret the CO2 spectrum include 
CO, each fit in part predicts a corresponding CO stretch- 
ing mode spectrum. This information is used to con- 
struct a model CO ice spectrum that can be compared 
directly with the observed 4.67 /im CO bands. The model 
CO profiles are therefore only fitted to the CO2 bending 
mode profiles and are plotted on the observed CO pro- 
files for comparison only. In summary, each 15.2 /im CO2 
bending mode is fitted with a function with 6 free param- 
eters; the depth of the 5 components and the mixing ratio 
of CO2 to CO for the blue component. The spectra and 
best fits are shown in Figures [3] to [THl along with the 
CO stretching mode bands, where available. The CO2 
column densities are given in Tabled 

It is found that the "blue" CO stretching mode band 
corresponding to the "blue" C02:C0 mixture was in 
general much too deep to fit the observed CO bands 
and it has consequently been reduced by a factor of 3 
relative to the CO2 bending mode in the comparison 
plots for all the sources. This discrepancy is discussed 
in Section 15.31 In this paper, the same band strength 

^ http://cow.physics.wisc.edu/~craigm/idl/idl.html 
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Fig. 3. — Sketch of the 5 different components used to fit the CO2 band. The spectra used to illustrate the components are from left to 
right: IRS 51, SVS 4-5 and RNO 91. 



of 1.1x10^^^ cm molecule^^ is used for ever y component 
of CO2. Note that while Gerakine s et alJ (JT995) mea- 
sure a larger band strength for the CO2 bending mode 
in a water-rich environment of 1.5xl0~^^cmmolecule~^, 
they also state that a large uncertainty is associated with 
this measurement. The effect of using this band strength 
for CO2 in a water-rich mixture would be to decrease the 
CO2 column density of the red component by 36%. 

5. RELATIONS OF THE CO2 COMPONENTS 

5.1. The abundance of CO2 ice in low-mass YSO 
envelopes 

The relation between the observed H2O ice column 
densities and the total CO2 ice column densities are 
shown in Figure [TH Because it is difficult to measure 
the column density of H2 gas along the line of sight, this 
is the relation typically used to define ice abundance of 
a given solid-state species as a number fraction relative 
to water ice. The CO2 ice abundance is remarkably con- 
stant, but does exhibit a scatter that is much larger than 
the uncertainties. A linear fit to the low-mass stars re- 
veals a number ratio of CO2 to H2O to be 0.32 ± 0.02, 
with a Pearson correlation coefficient of 96%. There is, 
however, a significant scatter in the relation, and a num- 
ber of our sources show abundances between 0.2 and 0.3, 
relative to H2O. The relation exhibits a slight tendency 
for higher CO2 abundances at higher H2O column den- 
sities. 

The CO2 abundance derived here can be compared 
with t hat of 0.17 for the ISO sample of massive 
YSOs (|Gerakines et all fl999l) and 0.18±0.04 for qui- 
escent clouds as o bserved toward background stars 
(|Wbittet et al. ' 2007). Both these samples have been in- 
cluded in Figure [ini The points associated with the mas- 
sive YSOs thus indicate a significantly lower CO2 abun- 
dance. Inspection of the figure also reveals that the back- 
ground stars generally probe lower column densities than 
the YSOs. For these low column densities, the difference 
between the YSOs and the background stars is less sig- 
nificant. However, at water ice column densities higher 
than ~ 2 X lO^^cm"^, the difference in CO2 abundance 
between background stars at low Ay and low-mass YSOs 
is highly significant. This sudden change in the relation 
may represent the activation of a new formation route to 
CO2 in addition to that forming the Il20:C02 mantle at 



lower Ay. An increase of CO2 abundance in den ser re- 
ions of a single cloud core was observed by Pontop pidanl 
..20061). and was a lso di scussed for the background stars 
in IWhittet et al.1 ()2007f ). based on a single star at high 
Av. 

5.2. The C02:H20 system 

The C02:H2 component dominates every CO2 bend- 
ing mode band, and the abundance variation of this com- 
ponent, or lack thereof, therefore mimics that of the to- 
tal CO2 band discussed in section [01 In the following, 
"abundance" refers to an observed column density ratio 
averaged over the line of sight, while "concentration" is 
the local point number density of a species relative to 
the total number density of molecules in an ice film. The 
abundance of this component varies between 0.1 and 0.3 
relative to water, with a median value around 0.2. If the 
concentration equals the abundance, the relevant mix- 
ing ratio for a laboratory analog is C02:H20=l:(5l2)- 
The observed abundances of the C02:H2 component 
are shown in Figure [20] as a function of the correspond- 
ing C0:H20 component of the CO stretching mode. 

The general value of this parameter is of some im- 
portance. [Oberg et al.l ()2007a[ ) showed that the band 
strengths of the various H2O modes are very sensitive to 
the CO2 concentration, and they suggested that this may 
be an explanation for the discrepancy in observed band 
depths between the H2O stretching and bending modes. 
The observed abundances of the C02:H20 component 
suggest that the influence of CO2 may explain a part of 
the H2O bending/stretching mode discrepancy. This is 
discussed in greater detail in Paper I. 



TABLE 4 
Ice column densities of the CO2 components" 



Source 


total CO2 


C02:H20 


C02;C0 ~1:1 


C02;CO~l:25 


Pure CO2 


shoulder 


C0:C02 ratio 


' X^ 


H20'= 


W3 IRS5 


6.56 ± 0.20 


4.41± 0.25 


0.93 ± 0.20 


0.02 ± 0.08 


0.83 ± 0.12 


0.27 ± 0.02 


0.82± 0.052 


0.49 


56.5 ± 6.0 


L1448 IRSl 


2.14 ± 0.06 


1.46± 0.11 


0.32 ± 0.10 


0.06 ± 0.03 


0.19 ± 0.06 


0.09 ± 0.01 


>1.00 


1.41 


4.7 ± 1.6 


L1448 NA 


40.92 ± 0.35 


32.76± 0.46 


3.79 ± 0.23 


0.00 ± 0.01 


4.21 ± 0.12 


0.43 ± 0.01 


<0.16 


1.82 




L1455 SMMl 


63.48 ± 4.43 


45.63± 1.78 


5.72 ± 1.68 


0.17 ± 1.06 


7.96 ± 0.49 


1.57 ± 0.05 


<0.16 


0.45 


182.0 ± 28.2 


RNO 15 


2.57 ± 0.05 


2.10± 0.08 


0.43 ± 0.04 


0.02 ± 0.02 


0.00 ± 0.05 


0.01 ± 0.05 


>1.00 


1.25 


6.9 ± 0.6 


IRAS 03254 


8.86 ± 0.10 


4.63± 0.18 


2.09 ± 0.12 


0.00 ± 0.01 


1.51 ± 0.09 


0.55 ± 0.01 


0.62± 0.030 


1.32 


40.5 ± 3.7 


IRAS 03271 


15.37 ± 0.09 


10.65± 0.16 


2.56 ± 0.11 


0.04 ± 0.07 


1.68 ± 0.06 


0.52 ± 0.01 


0.41± 0.012 


2.11 


76.9 ± 17.6 


B la 


20.85 ± 0.14 


14.23± 0.25 


4.25 ± 0.19 


0.31 ± 0.10 


1.66 ± 0.11 


0.37 ± 0.01 


0.65± 0.009 


1.31 


104.0 ± 23.0 


B Ic 


84.55 ± 15.70 


68.10± 3.72 


13.95 ± 2.02 


0.00 ± 1.37 


0.00 ± 0.10 


2.50 ± 0.11 


>1.00 


0.35 


296.0 ± 57.0 


IRAS 03439 


3.32 ± 0.06 


2.23± 0.11 


0.55 ± 0.09 


0.00 ± 0.03 


0.37 ± 0.06 


0.17 ± 0.01 


0.76± 0.082 


1.63 


10.1 ± 0.9 


IRAS 03445 


7.07 ± 0.09 


5.31± 0.17 


1.75 ± 0.12 


0.00 ± 0.06 


0.00 ± 0.07 


0.02 ± 0.07 


0.52± 0.018 


1.12 


22.6 ± 2.8 


L 1489 


16.20 ± 0.09 


11.40± 0.16 


1.80 ± 0.12 


0.00 ± 0.07 


2.66 ± 0.07 


0.04 ± 0.07 


0.44± 0.036 


1.95 


47.0 ± 2.8 


DG Tau B 


5.40 ± 0.06 


3.64± 0.11 


1.09 ± 0.09 


0.00 ± 0.01 


0.56 ± 0.06 


0.12 ± 0.07 


>1.00 


0.86 


26.3 ± 2.6 


GL989 


6.11 ± 0.07 


3.97± 0.19 


1.37 ± 0.13 


0.00 ± 0.04 


0.51 ± 0.09 


0.30 ± 0.02 


0.64± 0.013 


2.53 


23.2 ± 1.1 


HH46 


21.58 ± 0.11 


16.99± 0.20 


2.25 ± 0.12 


0.00 ± 0.09 


1.90 ± 0.06 


0.49 ± 0.01 


0.37± 0.016 


1.34 


77.9 ± 7.3 


CED 110 IRS4 


12.26 ± 0.12 


8.34± 0.22 


2.10 ± 0.17 


0.01 ± 0.09 


1.49 ± 0.10 


0.29 ± 0.02 


0.56± 0.009 


1.17 




B 35 


4.90 ± 0.15 


3.48± 0.27 


0.38 ± 0.20 


0.09 ± 0.09 


0.56 ± 0.14 


0.37 ± 0.02 


0.70± 0.223 


1.29 




CED 110 IRS6 


14.30 ± 0.08 


11.49± 0.14 


1.92 ± 0.08 


0.01 ± 0.07 


0.78 ± 0.04 


0.07 ± 0.04 


0.28± 0.005 


1.90 


47.0 ± 6.0 


IRAS 12553 


5.98 ± 0.09 


4.84± 0.16 


1.07 ± 0.11 


0.00 ± 0.05 


0.04 ± 0.07 


0.04 ± 0.07 


0.63± 0.014 


0.59 


29.8 ± 5.6 


ISO Chall 54 


1.81 ± 0.13 


0.50± 0.25 


0.96 ± 0.21 


0.17 ± 0.07 


0.16 ± 0.13 


0.22 ± 0.03 


>1.00 


1.23 




IRAS 13546 


8.72 ± 0.12 


6.22± 0.21 


2.00 ± 0.16 


0.01 ± 0.09 


0.30 ± 0.09 


0.08 ± 0.01 


0.54± 0.008 


0.86 


20.7 ± 2.0 


IRAS 15398 


52.16 ± 0.79 


38.83± 0.90 


0.29 ± 0.67 


1.85 ± 0.59 


10.58 ± 0.29 


0.80 ± 0.02 


0.16± 0.000 


1.34 


148.0 ± 39.5 


GSS 30 IRSl 


3.28 ± 0.06 


1.86± 0.10 


0.70 ± 0.07 


0.00 ± 0.05 


0.61 ± 0.05 


0.09 ± 0.05 


0.77± 0.078 


1.04 


15.3 ± 3.0 


WL 12 


4.34 ± 0.05 


2.72± 0.10 


1.36 ± 0.08 


0.04 ± 0.03 


0.03 ± 0.05 


0.20 ± 0.05 


>1.00 


1.28 


22.1 ± 3.0 


GY224 


2.69 ± 0.09 


1.90± 0.17 


0.66 ± 0.12 


0.00 ± 0.01 


0.03 ± 0.09 


0.15 ± 0.04 


0.65± 0.057 


1.28 




WL 20s 


5.02 ± 0.06 


3.86± 0.11 


0.75 ± 0.08 


0.00 ± 0.02 


0.28 ± 0.05 


0.11 ± 0.05 


0.64± 0.038 


0.95 




IRS 37 


4.05 ± 0.08 


2.99± 0.14 


0.99 ± 0.12 


0.00 ± 0.05 


0.01 ± 0.07 


0.01 ± 0.05 


0.65± 0.012 


0.85 


36.5 ± 5.0 


IRS 42 


4.49 ± 0.05 


3.39± 0.11 


1.01 ± 0.09 


0.07 ± 0.03 


0.01 ± 0.05 


0.09 ± 0.05 


>1.00 


1.01 


19.5 ± 2.0 


WL 6 


9.33 ± 0.08 


6.86± 0.15 


2.17 ± 0.10 


0.00 ± 0.06 


0.12 ± 0.06 


0.18 ± 0.05 


0.49± 0.007 


0.87 41.7 ± 6.0 


CRBR 2422.8-3423 


10.54 ± 0.06 


7.30± 0.11 


2.85 ± 0.07 


0.09 ± 0.04 


0.02 ± 0.04 


0.36 ± 0.05 


0.48± 0.003 


1.52 


45.0 ± 5.0 


IRS 43 


12.26 ± 0.12 


8.11± 0.23 


1.91 ± 0.16 


0.00 ± 0.09 


1.71 ± 0.10 


0.51 ± 0.01 


0.53± 0.047 


0.84 


31.5 ± 4.0 


IRS 44 


6.92 ± 0.08 


5.01± 0.14 


1.10 ± 0.11 


0.00 ± 0.01 


0.50 ± 0.08 


0.34 ± 0.01 


0.87± 0.040 


1.17 


34.0 ± 4.0 


Elias 32 


4.87 ± 0.09 


2.89± 0.15 


1.38 ± 0.12 


0.07 ± 0.05 


0.55 ± 0.07 


0.16 ± 0.01 


0.62± 0.015 


2.21 


17.9 ± 2.6 


IRS 46 


2.35 ± 0.12 


1.64± 0.23 


0.43 ± 0.20 


0.02 ± 0.07 


0.09 ± 0.12 


0.09 ± 0.08 


>1.00 


0.20 


12.8 ± 2.0 


VSSG 17 


5.86 ± 0.11 


3.46± 0.19 


2.20 ± 0.14 


0.00 ± 0.07 


0.01 ± 0.09 


0.27 ± 0.03 


0.53± 0.009 


0.64 


17.0 ± 2.5 


IRS 51 


9.32 ± 0.07 


5.44± 0.12 


3.30 ± 0.09 


0.27 ± 0.05 


0.07 ± 0.05 


0.26 ± 0.05 


0.54± 0.005 


0.88 


22.1 ± 3.0 


IRS 63 


6.84 ± 0.05 


4.49± 0.10 


2.26 ± 0.05 


0.01 ± 0.04 


0.00 ± 0.05 


0.17 ± 0.05 


0.51± 0.014 


1.20 


20.4 ± 3.0 


L 1689 IRS5 


3.37 ± 0.10 


2.19± 0.17 


1.05 ± 0.12 


0.00 ± 0.06 


0.00 ± 0.08 


0.22 ± 0.01 


0.58± 0.015 


1.26 




RNO 91 


11.66 ± 0.16 


5.94± 0.28 


2.37 ± 0.20 


0.00 ± 0.03 


2.70 ± 0.14 


0.62 ± 0.02 


0.64± 0.026 


1.49 


39.0 ± 5.0 


W33A 


14.11 ± 0.29 


9.96± 0.30 


1.70 ± 0.22 


0.02 ± 0.11 


1.13 ± 0.13 


1.41 ± 0.11 


0.59± 0.025 


1.00 


113.0 28.3 


GL 2136 


0.93 ± 0.03 


0.42± 0.03 


0.26 ± 0.02 


O.OOit 0.00 


0.15± 0.01 


0.16± 0.02 


0.57± 0.023 


0.69 


47.2± 4.7 


S68N 


43.27 ± 16.54 


30.99± 0.75 


5.38 ± 1.18 


0.09 ± 1.01 


6.81 ± 0.38 


0.00 ± 0.50 


<0.16 


1.90 




EC 74 


2.89 ± 0.08 


2.00± 0.14 


0.79 ± 0.07 


0.00 ± 0.05 


0.00 ± 0.01 


0.19 ± 0.01 


0.44± 0.032 


1.90 


10.7 ± 1.8 


SVS 4-5 


17.21 ± 0.10 


12.37± 0.18 


3.48 ± 0.13 


0.00 ± 0.08 


0.68 ± 0.07 


0.73 ± 0.02 


0.43± 0.009 


1.14 


56.5 ± 11.3 


EC 82 


2.54 ± 0.04 


1.88± 0.08 


0.07 ± 0.07 


0.12 ± 0.02 


0.34 ± 0.04 


0.08 ± 0.05 


>1.00 


2.18 


3.9 ± 0.7 


EC 90 


5.44 ± 0.05 


3.53± 0.09 


1.67 ± 0.07 


0.12 ± 0.03 


0.01 ± 0.05 


0.17 ± 0.02 


0.81± 0.011 


1.39 


16.9 ± 1.6 


SVS 4-10 


8.25 ± 0.05 


5.37± 0.09 


2.35 ± 0.07 


0.07 ± 0.04 


0.12 ± 0.04 


0.35 ± 0.01 


0.57± 0.004 


1.51 


16.0 ± 1.4 


CK4 


1.98 ± 0.09 


1.20± 0.16 


0.62 ± 0.08 


0.00 ± 0.05 


0.00 ± 0.05 


0.15 ± 0.00 


0.75± 0.099 


1.00 


15.6 ± 15.6 


CK 2 


11.93 ± 0.21 


9.02± 0.38 


2.29 ± 0.23 


0.01 ± 0.16 


0.26 ± 0.15 


0.31 ± 0.01 


0.32± 0.030 


0.75 


35.7 ± 3.5 


RCRA IRS5 


14.28 ± 0.13 


8.38± 0.23 


4.52 ± 0.17 


0.01 ± 0.10 


0.80 ± 0.10 


0.68 ± 0.02 


0.58± 0.005 


1.86 


37.6 ± 2.8 


RCRA IRS7A 


19.64 ± 0.12 


13.17± 0.21 


3.89 ± 0.16 


0.00 ± 0.09 


2.15 ± 0.08 


0.52 ± 0.01 


0.58± 0.006 


1.20 


109.0 ± 19.2 


RCRA IRS7B 


26.74 ± 0.22 


19.02± 0.35 


2.38 ± 0.24 


0.00 ± 0.25 


4.94 ± 0.11 


0.36 ± 0.10 


<0.16 


2.26 


110.0 ± 19.7 


IRAS 32 


18.70 ± 0.21 


10.21± 0.36 


4.00 ± 0.28 


0.13 ± 0.16 


3.04 ± 0.16 


1.32 ± 0.12 


0.56± 0.014 


1.52 


52.6 ± 18.8 


S 140 


3.78 ± 0.07 


1.17± 0.10 


1.06 ± 0.07 


0.00 ± 0.00 
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0.44 ± 0.13 
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NGC 7538 IRS 9 


15.81 ± 0.13 
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3.03 ± 0.12 


0.01 ± 0.06 


2.51 ± 0.08 
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0.89 


64.1 ± 6.4 
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32.51 ± 0.30 


23.66± 0.42 


4.20 ± 0.20 


0.00 ± 0.02 


4.10 ± 0.13 


0.40 ± 0.01 


<0.16 


3.02 


130.0 ± 22.6 



" All column densities are in lO^^cin^^. 

'' All uncertainties are statistical, and do not include systematic uncertainties from e.g. the continuum determination. 

'^ For consistency, the new water ice column densities from Paper I are used, where available. They are generally consistent with the 
few published values measured on the same, or similar, data sets, but with a few values differing by 15-20%. 
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Fig. 4. — Left panels: Continuum determinations for the SH spectra. The dashed curves are low-order polynomial fits, while the solid 
curves are the fits with a Gaussian included to simulate the blue wing of the 18 /xm silicate bending mode. Middle panels: Component fits 
of the CO2 bending mode profiles. Right panels: Th e CO stretching mode profiles predicted by the fit to the CO2 bending mode (where 
available) . The components are discussed in Section 14.31 The shaded areas indicate the profiles allowed by the uncertainties in the fit to 
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Fig. 18. — As Figure l4l but for archival ISO spectra of massive young stars from lGerakines et al.l II1999I) . 
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Fig. 19. — Relation between H2O ice column density and 
CO2 ice column density. The green points are low-mass stars 
with luminosities of 0. 1-100 L0, while the blue stars indicate the 
massive YSOs from the ISO sample with luminosities of L = 
10* — 10^ Lp, . The red t riangl es are the background s tars in Tau- 
rus fr om [Whittet et al.l | |200'7I') . as well as CK2 from IKnez et al.l 
ll2005l) .The dashed line is a linear least squares fit to the low-mass 
stars. The solid curve is the linear fit to the background stars. 
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Fig. 20. — Correlation between the C02;H20 component and the 
red CO component (C0:H20) from Pontoppidan et al. (2003b). 
The green points are the low-mass stars, while the blue points are 
the high mass YSOs observed with ISO. The line is a linear fit to 
the low-mass YSOs. 

5.3. The C02:C0 system 

5.3.1. The "blue" component 

The component fit indicates that the C02:C0 system 
forms a component separate from the CO2 mixed with 
water. The evidence is the higher abundance of the CO 
dominate d component in dens e, cold cloud regions as dis- 
cussed in iPontoppidanI ()2006[ ) , and as suggested by the 
offset of low-mass protostellar envelopes relative to back- 



ground stars and massive YSOs in Figure [TH as well as 
the match of the profile of this blue CO2 component with 
laboratory simulations of C02:C0 mixtures. An impor- 
tant parameter to determine is the concentration of CO2 
relative to CO in this component. There are at least 
two ways of doing this. First, the concentration can be 
directly determined by comparing the observed strength 
of the blue CO2 component with the corresponding blue 
component of CO observed as part of the CO stretching 
mode band. Second, it can be indirectly inferred by the 
profile of the C02:C0 component of the CO2 bending 
mode, since this is sensitive to the concentration. 

The column density ratios of the CO2 and CO 
"blue" are illustrated in Figure [5T] These compo- 
nents exhibit a fairly strong correlation with a Pear- 
son correlation coefhcient of 0.70 and a slope of 
iV(blue C02)/iV(blue CO) = 2.5 ± 0.2, assuming a width 
of 3cm~^ f or the blue CO component as empirically de- 
termined in iPontoppidan et al.l (|2003bt ). The laboratory 
spectra of the CO stretching mode of C02:C0 mixtures 
are about twice as wide. Using the laboratory spectra 
instead to calculate the column densities of the "blue" 
CO component would decrease this ratio to 1.25±0.1. 

Conversely, Figure [H] illustrates the indirect method 
of determining the concentration of CO2 in CO from the 
profile of the "blue" CO2 component. The figure shows 
the distribution of C02:C0 mixing ratios as determined 
by the component fit. It is seen that the mixing ratios 
are remarkably similar for most of the observed spec- 
tra, which is also indicated by the direct correlation seen 
in Figure [21] However, the median concentration deter- 
mined using this method is iV(blueC02)/Af(blueC0) = 
0.55±0.2, or a factor of 2-5 smaller than that determined 
using the directly measured column densities. 

It is probably reasonable to assume that the direct 
method provides a better estimate of the concentration 
since the indirect method relies on an uncertain calibra- 
tion of a set of laboratory experiments. However, it 
should be noted that the band strengths of both CO2 
and CO may depend on concentration, which is an effect 
that is ignored here by assuming that the band strengths 
are constant. Variable band strengths may affect the di- 
rect method. Clearly, well-calibrated laboratory exper- 
iments are needed to resolve the issue. For this study, 
the "blue" CO components, as determined by the profile 
of the blue CO2 band, are divided by a factor 3 to bet- 
ter match the CO bands, as dictate d by the direct con- 
centrat ion measurement. Note that IPontoppidan et all 
(|2003bf) in their study of the CO stretching vibration 
band at 4.67 /xm found that the available laboratory spec- 
tra of C02:C0 mixtures were generally much too wide 
to reprod uce the blue wing of the CO bands, as con- 
firmed bv van Broekhuizen et al.] (|2006f ). which led them 
to consider alternatives to explain the "blue" component 
of the CO band. Here it is concluded, based on the clear 
correlation between the blue CO2 and CO components, 
that they indeed represent the same component, but that 
since their detailed profiles do not fit well to the labora- 
tory analogs, the structure of this component is not yet 
fully understood. 

5.3.2. The "dilute" component 

Some CO2 bending mode profiles exhibit a very narrow 
single peak at 15.15 /im (660 cm~^), as indicated in Fig- 
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Fig. 21. — Correlation between the blue CO2 compo- 
nent abundance and the blue CO component abundance from 
IPontoppidan et al.l l l2003bl ). otherwise as Figure [20l 
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Fig. 22. — Distribution of C02:C0 mixing ratios as determined 
from the CO2 blue component profile fit. 



ure [31 The most obvious sources with this property are 
IRS 51 and CRBR 2422.8-3423, as well as the background 
star CK 2. The profile of this component corresponds 
closely to that of a CO2 ice very dilute in CO with a con- 
centration less than 10%; the profile of the band does not 
change appreciably at lower concentrations. The "dilute" 
component typically appears toward sources that also 
have very large col umn densities of "pu r e" CO (the "mid- 
dle" component of IPontoppidan et al.1 (|2003b[ )). The re- 
lation between these two components is shown in Figure 
[221 where it is seen that typical concentrations of CO2 in 
the CO is 1:100-250. This indicates that there are vastly 
different mixing ratios of CO2 to CO along each line of 
sight, possibly even on each grain. 



0.020 




0.000 



0.5 1.0 

N(Pure CO)/N(H^O) 

Fig. 23. — Relation between the dilute CO2 component abun- 
dance and the middle Jjpuve) CO component abundance from 
Pontoppidan et al. (20032)- The green points are detections, while 
the red points indicate upper limits on the dilute CO2 component. 
The line is a linear fit to the detections. 

5.4. Relation with CH3OH 

Because CH3OH has been related to th e shoulder on 
the r ed side of the CO2 bending mode (jDartois et al.l 
|1999( ). it is natural to estimate whether there is a con- 
nection with the direct measurement of CH3OH abun- 
dances from Paper I. The relation, shown in Figure [Ml 
does not exhibit an obvious relationship between the CO2 
shoulder and the CH3OH abundance. This does not nec- 
essarily indicate that the shoulder is not related to in- 
teractions with CH3OH if the concentration of CO2 in 
the CH3OH varies significantly. Keeping the assumption 
that the band strength of the CO2 shoulder is that of 
pure CO2, the abundance C02:CH30H varies between 
1:20 and 1:3. It is therefore likely that the CO2 is highly 
dilute in the CH3OH ice. 

5.5. Upper limits on C3 O2 

iJamieson et al.l (|2006[ ) found that carbon suboxide 
(C3O2) is formed in abundance along with CO2 during 
electron irradiation of a CO ice. I t has prominent bands 
around 17-19 /im (Ger akines fc M oore 2GGla), but with 
exact central positions that vary considerably in the lit- 
erature, presumably as a result of different ice mixtures. 
Some of the IRS spectra do have weak structure in the 
general wavelength region, but nothing that resembles a 
consistently recurring absorption band at a single wave- 
length. It is therefore concluded that there is no obvious 
evidence for absorption from C3O2 at the 5% level. 

5.6. The 17 ^m feature 

Some spectra show a clear excess absorption feature 
centered on 17/J,m, most notably toward IRS 42, EC9G, 
VSSG 17 and ISO Chall 54. The origin of this feature 
is unclear, and only its observed properties are reported 
here. It appears that there is some relation between the 
presence of the feature and sources that show contami- 
nation by silicate emission from the central disk. Thus 
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Fig. 24. — Relation between the abundance of the CO2 "shoul- 
der" component and the CH3OH abundance from Paper I. 

the feature may not be a real absorption signature, but 
emission from the 18 /im sihcate band affecting the con- 
tinuum determination. However, if real, the feature has 
a typical width of 35cm~^, and optical depths of 0.1, 
where detected. 

6. ORIGIN AND EVOLUTION OF THE CO2 COMPONENTS 

6.1. Formation routes to CO2 in the interstellar 
medium 

CO2 ice in protostellar envelopes as observed here 
is widely believed to be formed by surface reactions 
as opposed to a gas- phase route followed by freeze-out 
(jTielens &: Hageij TT982V Pure gas-phase chemical mod- 
els of typical cold, dark clouds predict CO2 abundances 
of 10~^ relative to H2 (Bergin ct al. 1995), making it un- 
likely that there is any contribution from directly from 
gas-phase chemistry to the observed sohd state CO2. 
There is strong observational evidence that CO2 forms 
readily in cold quiescent molecular clouds, and it does 
not appear that "extra" photo-processing of the ice is 
required, beyond what can be explained by cosmic ray 
induced UV photons and photons from the interstel- 
lar radiation fi eld managing to pen etrate to Ayn of 3- 
5 magnitudes (jWhittet et al.|[T998[ ). In fact, the rela- 
tively constant fraction of CO2 relative to water ice (15- 
40%) under a wide range of conditions suggests that it 
forms under "common" quiescent conditions of densities 
of 10^ — 10^ cm^'^, temperatures near 10 K and a standard 
cosmic ray field. Thus far, no line of sight has been ob- 
served to have CO2 abundances of A^c02/-^H20 ^ 0.15. 
However, the exact chemical pathway remains controver- 
sial. Possible routes to the formation of CO2 are via the 
reactions: 

CO + OH^COa+H (1) 



O -H HCO -* CO2 + H 
C0-K0^C02 



(2) 
(3) 



120041) ■ but has also been found in at least one 
study to possess a pro hibitively large barrier 
(|Grim fc d'HendecourtI 119861). However, a similar 
experiment bv lRoser et al.l (|2001h finds that the reaction 
proceeds with no or little barrier. It should be noted 
that it is expected that the barrier to route [3] is sensitive 
to the electronic state of the oxygen atoms, such that 
0{}S) may react much ea sier with CO than oxygen 
in the ground state 0{^P) (jFournier et al.lll979f ). The 
energy difference between these two states correspond 
to a red photon (6300 A), which will penetrate much 
deeper into dark clouds than the UV photons normally 
considered for photolysis reactions. 

In some grai n surface models, ro ute[T]is used as a dom- 
inant reaction (jChang et al.ll2007l ). 

It is also known that CO2 can form with a low or non- 
existing activation barrier through an electronically ex- 
cited state of CO: 



CO* + CO ^ CO2 + C 



(4) 



Route 13 is often included as an i r nportant grain surface 
reaction (jTielens fc Hageiil 119821 : IStantcheva fc Herbstj 



Obviously, this reaction requires that the CO molecule 
is excited, and was studied extensively in the context of 
UV photolysis (e.g. iGerakines et"anil996l : lLoeffler et al.l 
l2005f ). However. I Oberg et al.l (|2007bD finds no formation 
of CO2 from CO in a ultra-high vacuum (UHV) UV irra- 
diation experiment, which is much less contaminated by 
H2O than the previous high vacuum experiments. While 
this argues against route[4]as an effective pathway to CO2 
ice, electron irradiation does produce CO2 from pure CO 
under UHV conditions (Jamieson et al. 2006) . 

Consequently, the rates of most solid-state reactions 
leading to CO2 are still controversial, and theoretical 
models have struggled to consistently reproduce the ob- 
served abundances. Based on existing observations, any 
model for the formation of CO2 should be required to re- 
produce both the absolute abundance of CO2 of 15-40% 
relative to water ice, as well as the apparent universality 
of this abundance. Additionally, the separate molecular 
environments of the CO2 ice should also be explained, 
in particular the presence of CO2 in both H20-rich and 
CO-rich environments. 

It is worth mentioning that CO accreted from the gas- 
phase is not the only p otent ial source for the carbon in 
CO2. iM ennella et al.l (|2004| ) found that the carbon in 
hydrogenated carbon grains will form both CO and CO2 
when covered with a water ice mantle and subjected 
to cosmic rays. This is potentially a route to forming 
CO2 embedded in a water ice matrix. However, hav- 
ing a different source of carbon for CO2 formation than 
gas-phase CO must explain why the ^^C/^'^C ratios of 
gas-phase CO, solid state CO and solid st ate CO2 are all 
so similar to the Solar value of 89 (50 -lOO lBoogert et al.l 
I2OOOII2OO2I : iPontoppidan et al.ll2003"bl ). In contrast, pre- 
solar carbonaceous grains have highly variable ^^C/^'^C 
ratios with a tendency toward ratios a few to 100 times 
that of the Sun, although some have ra tios as low as ~1 
(e.g. lLin et all 120021 : ICroat et al.ll2005( ). The question is 
whether such scatter is reflected in the ice isotopologue 
ratios if the carbonaceous grains are the source of carbon. 

In the following, it is explored how the data presented 
here can help constrain the formation and evolution of 
CO2 in the interstellar medium and in protostellar en- 
velopes in particular. 
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6.2. Formation of the C02'-C0 system 

Having established its existence, how can the presence 
of large amounts of CO2 within the CO-doniinated man- 
tle be explained? Is there a connection with the freeze- 
out of CO at densities higher than 10^ cm^'^? There is 
evidence from laboratory experiments that CO and CO2 
do not mix upon warmup of separately deposited layers 
(|van Broekhuizen et al.l [20061 ). and the possibility that 
the CO2 is formed directly as part of the CO mantle us- 
ing the carbon of the CO accreted from the gas-phase 
is therefore explored. While a thick mantle of CO ice 
certainly offers a significant reservoir of carbon and oxy- 
gen for the formation of CO2, the breaking of the CO 
triple bond to form CO2 directly from CO requires a 
significant energy input. The most direct way of pro- 
viding a high input of energy in dense molecular clouds 
is via cosmic rays. The cosmic rays can hit the grains 
directly, but are also the dominant source of UV pho- 
tons through their interaction with hydrogen molecules. 
A number of laboratory experiments have been per- 
formed simulating and comparing the formation of CO2 
from a pure CO ice through UV and cosmic ray irra- 
diation ilGe rakines fc Moore"20Glb'; 'Loefflcr e t al1l2005l : 
iJamiesoiiet al. 2006). Ja micson et al. (20061) showed 
that CO2 can be formed from a pure CO ice layer during 
irradiation with a 5 keV electron beam simulating heav- 
ier and more energetic cosmic rays. The experiment con- 
verted 0.49% of the CO to CO2 with a deposited energy 
of 5.8 X 10^^ MeVcm"'^, or 3.4eVmolecule~^, assuming 
a CO ice density of 0.8gcm~'^. 

This value can be put into the context of dense molec- 
ular clouds by estimating the time scale for depositing 
the same amount of energy to a CO ice mantle with 
a standard cosmic ray field. Following the approach of 
IShen et all (|2004f ). the total deposited energy by the cos- 
mic ray field per CO molecule per second is: 



Edc-piCR) 



mco v^ 



Pco 






47r 



,{Z) as dE 



(5) 



where toco is the mass of a CO molecule and pco is the 
density of CO ice. dQ/ds is the energy loss as the cos- 
mic ray traverses a length, ds, through the CO ice and 
dn/dE is the cosmic ray flux spectrum, and finally, fz is 
the fraction of the energy deposited that actually remains 
with the g rain the rest being ejected in highly energetic 
electrons. iLeger et"al] ()1985D estimate that fz ~ 0.6. 
The minimum cosmic ray energy, i^min, is determined by 
the various factors that may drain energy from the par- 
ticles. These include interactions with dust grains, with 
the molecular gas itself, as well as drag from magnetic 
fields. The importance of these effects depend sensitively 
on Z as well as the cosmic ray energy. iLeger et"al] ()1985f ) 
showed that interactions with the gas only affects the cos- 
mic ray spectrum at Ay > 50 mag. However, dust can 
effectively stop low energy (< 100 MeV for iron) cosmic 
rays at Ay < 10 mag, and that determines i^min- 

For these assumptions, the energy deposited by cosmic 
rays in a CO ice mantle is 8 x 1 ~^^ eV molecule"^ s~'^, 
which means the iJamieson et al.l ()2006[ ) experiment cor- 
responds to roughly 350 yrs of cosmic ray irradiation or, 
assuming a constant rate of CO2 formation, ~ 3.5 x 
10** yrs to convert 50% of the CO mantle to CO2. Al- 



though the input cosmic ray flux spectrum is very uncer- 
tain, especially at low energies, the conclusion is that it 
is plausible that cosmic rays can provide the necessary 
energy input to form the observed C02:C0 component. 

In this context, what is the implication of the presence 
of the "dilute" component discussed in Section I5.3.2P 
One possibility is that the CO ice with a low concen- 
tration of CO2 is younger than the C02:C0~1:1 compo- 
nent. This would happen if the conversion of CO to CO2 
o ccurs at a con stant rate. 

IJamieson et a l. (2006) predicts the presence of a range 
of carbon oxide species in addition to CO2, of which the 
most abundant is C3O2. The quality of the Spitzer-IRS 
spectra allows a sensitive search for this molecule in the 
solid state through its modes at ~ 18.5 /zm, yet it is not 
clearly detected in any of the spectra presented here (see 
Section [53]). 

6.3. CO2 as a temperature tracer 

CO2 has been suggested to be a tracer of strong heat- 
ing based on simula ted annealing experim ents in the 
laboratory to 100 K ([Gerakines et al.[[1999D . The pro- 
posed mechanism is that the CO2 segregates out of 
the hydrogen-bonding mixture with water and possibly 
CH3OH to produce inclusions of pure CO2. These in- 
clusions in turn produce the characteristic double peak 
observed in many high mass YSOs. While the very 
high temperatures required for the segregation process 
in a laboratory setting probably correspond to somewhat 
lower temperatures on astronomical time scales, they are 
still well above the temperatures of 10-40 K that domi- 
nate the column densities through protostellar envelopes 
around low-mass stars. Thus, it seems surprising that 
many of the surveyed low-mass stars show a double peak. 
The decomposition and visual inspection of the spectra 
reveals that a pure CO2 component is clearly detected in 
18 of the 48 low-mass stars, or almost 40%. 

At least one other mechanism to produce a pure CO2 
ice component exists. The existence of a ubiquitous 
C02:C0 component has been suggested before and is 
strengthened by the sample presented here. This com- 
ponent may also produce pure CO2 through distillation. 
Upon warmup of a C02:C0 mixture the CO will desorb 
leaving the CO2 behin d, but will do so at rnuch l ower 
temperatures (20-30 K van Broekhuizen et al.[[2006[ ). 

One way of distinguishing this process with the forma- 
tion of pure CO2 via segregation is to calculate the frac- 
tion, P, of thermally processed icy material at tempera- 
tures above a certain critical temperature, Tcrit, along 
a given line of sight through a protostellar envelope. 
Clearly, Tcrit is expected to be much lower for the distilla- 
tion process than for the segregation process. Assuming 
r is a monotonically decreasing function of radius, P is 
given by: 



P 



Ir'"^ rico2ir)dr 



rR.u 



loJ"" nc02ir)dr 



(6) 



where i?crit is the radius where T = T^it and i?sub is the 
radius where the processed CO2 ice sublimates. nco2 is 
the density of the CO2 component that is transformed to 
pure, crystalline CO2 upon heating to Tcrit. The use of 
a critical temperature assumes that the process forming 
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Fig. 25. — Relation between the CO and CO2 temperature trac- 
ers. The CO2 temperature tracer is the ratio of the column density 
in the double peaked, pure CO2 component relative to CO2 in wa- 
ter. The CO temperature tracer is the amount of pure CO relative 
to CO trapped in water ice. Green points are low-mass YSOs, blue 
points are high-mass YSOs, and red points are upper limits on the 
CO2 temperature tracer. The line is a fit to the low-mass YSOs. 

pure CO2 is a thermal process governed by some activa- 
tion energy and described by an Arrhenius relation. 

The first step is to estimate the value of Tcrit for the 
segregation process in water-rich ice. Because of the long 
time scales and low pressures in the interstellar medium 
compared to the short time time scales and high pres- 
sures of a laboratory experiment, it is not appropriate 
to apply the critical temperatures from the laboratory 
directly to an astrophysical problem. For a process not 
dependent on pressure, such as the segregation of CO2 
from a water matrix, the critical temperatures in the two 
settings are related via: 



'^astro 
Tlab 



exp 



E,. 
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T«, 



1 

Tiab 



(7) 



where Ea is the activation energy in Kelvin, while Tastro 
and Tiab are the e-folding time scales of a given process 
in the interstellar medium and in the laboratory, respec- 
tively. Because Ea is unknown, measurements of the 
laboratory time scale at two different temperatures are 
required. 

Unfortunately, the kinetics of the segregation pro- 
cess are not well known. From the experiment of 
lEhrenfreund et all (jl999j , a rough estimate can be made 
of Tiab ~ Ihr at 100 K and Tiab ~ Imin at 120 K. 
However, they use a tertiary C02:H20:CH30H=1:1:1 
mixture, which has a concentration of methanol much 
higher than that found in typical low-mass protostcUar 
envelopes. Conversely, recent experiments with a binary 
C02:H20=1:4 mixture by lOberg et al] (|200l3) show 
that the CO2 bending mode double peak has formed al- 
ready at 75 K on laboratory time scales of hours. Using 
a time scale of 10^ yr, a typical time for the collapse 
front to reach the outer boundary of the protostellar 
core, the Ehrenfrcund et al. values give E g ~ 4 900 K 
and Tcrit = Tastro ~ 70 K. iBoogert et all (|2000f) find 



Ea = 4900 and Tcrit = 77K with similar assumptions. 
If it is instead assumed that Tiab — 75 K for an e- 
folding time scale of 1 hour, as indicated bv lOberg et al.l 
(|2007al) . but the activation energy of 4900 K is retained, 
Tcrit ^ 57 K. It is stressed that these values for Tcrit are 
educated guesses at best, and that quantitative kinetic 
laboratory experiments are needed to measure the ac- 
tual value. In conclusion, Tcrit is taken in the range 60 
to 80 K. 

It is also important to note that the time a dust grain 
can be expected to spend at temperatures between, say, 
70 and 90 K is much less than 10^ years. In the sim- 
plest physical 1 -dimensional model of an infalling enve- 
lope (|Shulll977l ). a dust grain at the radii correspond- 
ing to such temperatures will be in free fall. The time 
scale for it passing through this region for a typical 1 Mq 
young star is 75 years, which in turn will increase Tcrit 
for segregation to 78 K for the CHsOH-rich mixture. Us- 
ing a 2-dimensional infall model that takes rotation into 
account will likely increase the infall time scale some- 
what. The confidence of the value of Tcrit can obviously 
be improved significantly with a quantitative laboratory 
simulation coupled with a more detailed infall model. 

A value for the desorption temperature, Tgub, of the 
pure, crystalline CO2 component formed by the segre- 
gation process is also needed. It is reasonable to expect 
the segregated CO2 to be in the form of inclusions em- 
bedded in the water ice. The question is whether the 
CO2 is trapped in the water or will be able to escape at 
temperatures lower than the interstellar water ice des- 
orption temperature of 110 K (Frascr et al. 2001). While 
the Ehrenfreund experiments retain the CO2 inclusions 
until the water ice desorbs at 150 K, the Oberg et al. 
experiments find that bulk of the CO2 ice desorbs at 
temperatures much lower t han the water i ce. T his is 
consistent with the result of ICoUings et al.l (|2004l ). who 
classified CO2 as a molecule that is not easily trapped 
in a water ice matrix. The fact that the tertiary mix- 
ture Ehrenfreund experiment appears to retain the CO2 
to higher temperatures than the binary II2O mixtures 
may be related to the CH3OII changing the trapping 
properties of the matrix, along with In any case, to 
match the different experiments, an interstellar CO2 des- 
orption temperature of 110 K is assumed for the Ehren- 
freund experiment, 60 K for the Oberg experiment and 
50 K for the pure CO2 layer produced by the distillation 
process of C02:C0. The physical interpretation is that 
the methanol-rich ice traps the CO2 ice until the water 
desorbs. The water-rich binary mixture only traps CO2 
to temperatures slightly higher than the desorptio n tem- 
peratu re of pure CO2 in accordance with CoUings et al.l 
(|2004f ). while C02:C0 mixture obviously does not trap 
CO2 at all. 

The next step is to choose a radial temperature and 
density structure that can be used to calculate P. While 
P is sensitive to a range of structural parameters for 
the envelope, the dominant one is the luminosity of the 
central source. It is beyond the scope of this paper to 
explore the parameter space of the structures of proto- 
stellar envelopes, but it is instructive to construct an 
example. For simplicity a static, one-dimensional power- 
law envelope (p(i?) esc (i?/i?o)~^'^) with Ay of 50 mag 
and i?o = 100 — 300 AU is assumed. Furthermore, it 
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is assumed that the envelope is empty within Rq. This 
is actually an envelope structure that favors a large P. 
More evolved envelopes dominated by infalling material 
to large radii will have a shallower density profile and 
thus more of the line of sight column density at larger 
radii. In this sense, the model P curves are upper lim- 
its. The dust temper ature is calculated using the M onte 
Carlo code RADMC iDullemond fc DominikI ()2004[ ) cou- 
pled with a dust opacity constructed to fit the extinction 
curve as measured using the c2d photometric catalogues 
(Pontoppidan et al., in prep.). The resulting dust tem- 
perature at 100 AU varies from 40 to 1200 K for source 
luminosities of 0.1 to 10^ L©. 

Figure [^Hl shows the observed values of P as a func- 
tion of source luminosity compared to model curves for 
different values of Tdit under the assumptions that the 
pure CO2 originates either in the C02:H2 component 
or in the C0:C02 component through the processes dis- 
cussed above. The P-curves for both the Ehrenfreund 
et al. and Oberg et al. experiments are shown, as well 
as the curves expected for the formation of pure CO2 
through distillation of the C02:C0 component. 

First, it is noted that P is not necessarily a mono- 
tonic increasing function with luminosity. This is be- 
cause the evaporation of the pure CO2 ice component in 
the innermost regions of the envelope where T > Tguh 
competes with the formation of pure CO2 at tempera- 
tures Tcrit < T < Tsuh- This is seen in Figure US] as a 
turnover in the models as the luminosity increases. For 
higher power law indices of the envelope, the P-curves 
may even decrease with increasing luminosity. 

Comparing the data points with the model curves, the 
results are as follows: Assuming a complete transforma- 
tion of CO2 mixed with water ice to pure CO2 inclusions, 
the observed points are consistent with a critical temper- 
ature for this process of 50-70 K, depending on whether 
the Ehrenfreund et al. or the Oberg et al. experiments 
are considered. Conversely, assuming a conversion of 
CO2 mixed with CO to pure CO2, a critical temperature 
of at most 25 K explains the highest observed P values 
toward low-mass stars. The measured values of P are 
very sensitive to the presence of cold foreground clouds 
contributing an unrelated column density, as well as the 
detailed structure of the inner envelope, in particular the 
arbitrary location of an inner edge at 100-300 AU. This 
has important consequences for the use of the the CO2 
bending mode as an astrophysical tracer. For instance, 
the models show that the double peak should have a 
roughly constant relative strength for sources with lumi- 
nosities between a few and at least 10'^ L©. Therefore, if 
a source within that luminosity range shows no sign of a 
double peak, it is an indication of the presence of a sig- 
nificant contribution to the extinction from foreground 
material, unrelated to the protostellar envelope. 

In conclusion, the splitting of the CO2 bending mode 
toward low-mass protostars can be explained by segre- 
ga tion in strongly heate d water-rich ices as described 
in lGerakines et aLl ()1999f ) only for protostellar envelopes 
with steep density density profiles extending all the way 
to 100 AU from the central star, and for concentrations 
of CII3OII much higher than the observed abundances. 
The experiments of lOberg et al.l (|2007af l may also explain 
the data for such envelopes, but Tcrit is not well-defined 



for this experiment. 

The resulting distribution of CO2 ice in the different 
environments is sketched in Figure [571 

7. DISCUSSION AND CONCLUSION 

In this paper, a picture is presented in which the dom- 
inant ice components, II2O, CO2, CO and, in a limited 
number of cases, CH3OII appear to constitute an inti- 
mately connected system; the observational characteris- 
tics of each species are directly affected by the presence 
of the others. Other molecules act as trace constituents 
with abundances that are too low to significantly affect 
the profiles of the dominant species, i.e. other species 
did not have to be invoked to explain the observed band 
profiles. This is in contrast to the very complex 5-8 /im 
region discussed in Paper I. 

The survey has shown that the CO2 bending mode pro- 
files toward low-luminosity (0.1 < L < 100 L©), solar- 
type YSOs do not differ strongly in their basic struc- 
ture from those observed toward massive, luminous YSOs 
(L > 1000 L©). All protostellar envelopes as well as qui- 
escent molecular clouds probed by lines of sight toward 
background stars are dominated by CO2 mixed with wa- 
ter with abundances of 20% relative to water ice, but 
with a significant additional contribution of 10% of CO2 
mixed with CO. The other components needed to explain 
the CO2 bending mode are present with much smaller 
abundances, or a few % each with respect to water ice. 

While the C02:Il2 component might form from CO 
and OH via reaction ([1]), it is probable that, because of 
the lack of residual atomic hydrogen, the C02:C0 com- 
ponent forms through a different chemical route. Addi- 
tionally, since the band profiles indicate that the CO2 
in this component can have concentrations of anything 
from 2:1 to less than 1:100 relative to CO, a variable 
effect highly sensitive to either the environment or age 
of the ice must play a role. Based on ob servations of a 
sample of background stars. iBergin et all (2005.) suggest 
that the formation of the C02:C0 component is related 
to the abundance of atomic oxygen relative to atomic 
hydrogen in the gas-phase and predict that it forms in 
low-density regions of the cloud. Observationally, this 
would produce CO and CO2 profiles dominated by the 
C02:C0~1:1 component in low-density regions. It would 
also result in the C02:C0 layer being placed below the 
water-rich layer on each grain. In this case, the forma- 
tion of pure CO2 by the desorption of CO as proposed 
here would not work. 

A diflFerent scenario is suggested in which the C02:C0 
component is connected to the rapid freeze-out of pure 
CO at higher densities. This scenario is practically inde- 
pendent on the gas-phase chemistry, but requires a mech- 
anism for forming CO2 directly from CO. Such a mech- 
anism will most likely involve a highly energetic input 
from cosmic rays; UV irradiation with A > 1200 A will 
not work in the absence of water (jOberg et al.ll2007H) . 
While the energy input from a standard cosmic ray field 
is likely sufficient, this process would also tend to form 
more complex carbon oxides in the ice. While these are 
are not detected in our Spitzer spectra, their absence 
is not strongly constraining, given that the molecular 
properties of their bands, such as strengths, arc not well 
known. 

It is confirmed that the CO2 ice profile is an excellent 
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Fig. 26. — Fraction of pure CO2 compared to models for different critical temperatures relevant for segregation (left) and CO desorption 
(right). In the models on the left and middle panels, an inner edge at 100 AU is used, while the models on the right panel fit better with 
an inn e r edge at 30 AU. The dots are the obs er ved values fo r the subse t of our sample t h at ha v e measured lum inosities from Bcrrilli et id] 
flpgl); 'Ladd ct al.' (Togf); 'Che n et al] ||T995| 'I: ISaraceno et al., (,19961) : IBontemps et al.l 1I2OOII) : ILarsson et"all (i200a) : .Kaas et al., (,20041) 7 

lfePteS'l^Sftife:P6''i?S6€¥^. 'SSmparison with the stretching create pure CO2 by distillation rather than segregation. 



^^m?j0S¥i\h^Pi^^!?SSS^. (ByMparison with the stretching 
vibration band of sohd CO shows that the CO2 corre- 
lates weU with the pure CO to C0:H20 ratio, an estab- 
hshed ice temperature tracer. The observed differences 
between CO2 bending mode profiles are consistent with 
differences in the fraction of the sight lines that have been 
heated above a certain threshold temperature, Tdit, to 
produce a double-peaked profile. More luminous YSOs 
have generally had a larger fraction of their ice column 
densities at temperatures above Tcrit, but the prevalence 
of the characteristic double-peak shows that even the 
low-luminosity YSOs have thermally processed inner en- 
velopes. 

There are several possibilities for the value of Tcrit- For 
the massive YSOs, strong annealing of methanol-rich ices 
to temperatures higher than 100 K in a laboratory setting 
was invoked _to explain the observed abundance of pure 
CO2 ice by Gerakines et al. (1999). It is suggested that 
this corresponds to Tcrit = 50 — 80 K under the conditions 
in collapsing protostellar envelopes. These are high tem- 
peratures that puts restrictions on the density structure 
of the protostellar envelopes in order to reproduce the 
observed abundances of pure CO2. Laboratory experi- 
ments measuring the kinetics of the annealing process in 
methanol-poor C02:H2 ice mixtures will be needed to 
use this process in quantitative modeling of ice mantle 
processing. 

It is therefore argued that the presence of a significant 
fraction of pure CO2 toward so many of the low-mass 
YSOs in the survey indicates that another, lower temper- 
ature process may play a role. Our survey, as well as oth- 
ers, measure a significant fraction of the CO2 ice mixed 
with CO, rather than with water. This mantle compo- 
nent will form a layer of pure CO2 upon very moderate 
heating to Tcrit =20-30 K as the CO desorbs, leaving the 
CO2 behind. This process has therefore the potential to 



Note that the new observations presented here do not 
rule out that the segregation process is responsible for a 
significant fraction of the pure C O2 in the sample of mas - 
sive YSOs originally discussed in lGerakines et al] ()1999f) . 
in particular since less CO will be frozen out in the warm 
envelopes surrounding massive YSOs. The simple proto- 
stellar models presented here seem to indicate that both 
mechanisms for producing CO2 in general contribute to 
the double-peak in low-mass young stellar objects. 
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